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1 Summary 

An  adaptive  control  methodology,  merging  two  known 
approaches  to  flight  control  problem,  gain-scheduling  and 
direct  eigenspace  assignment  (DEA),  is  developed.  A gain- 
scheduled  inner  (stability)  loop  structure  is  shown  to 
minimize  the  variance  of  the  outer  (guidance)  loop  gains  and 
increase  the  robustness  of  the  system.  The  employment  of 
DEA  with  gain  scheduling  is  observed  to  decouple  the 
longitudinal  and  lateral  flight  modes  resulting  adequate 
system  stability,  enhanced  robustness  and  control  surface 
effectiveness.  This  methodology  is  used  in  the  flight  control 
law  design  of  a UAV. 

2 Introduction 

Employment  of  eigenvalue  / eigenstructure  assignment 
techniques  has  found  widespread  practical  applications  that 
solve  a variety  of  flight  control  problems.  Digital 
implementation  of  these  applications  has  evolved  to  be 
commonplace  by  the  recent  advances  in  microprocessor 
technology. 

In  this  paper,  two  known  approaches  to  the  flight  control 
design,  gain  scheduling  and  eigenspace  assignment  are 
reviewed.  In  the  third  and  fourth  sections,  two  methods  are 
explained  with  their  formulations,  respectively.  In  section 
five  the  flight  control  problem  is  stated  and  developed 
scheme  is  presented  to  solve  the  problem,  to  demonstrate 
synergetic  use  of  these  two  methods.  Section  six  discusses 
the  results  achieved. 

3 Gain  Scheduling 

In  many  situations  it  is  known  how  the  dynamics  of  a process 
change  with  the  operating  conditions  of  the  process.  One 
source  for  the  change  in  dynamics  may  be  the  nonlinearities 
that  are  known.  It  is  then  possible  to  change  parameteres  of 
the  controller  by  monitoring  the  operating  conditions  of  the 


process.  This  idea  is  called  gain  scheduling.  It  can  be 
regarded  as  a special  kind  of  nonlinear  open-loop 
adaptation  of  regulator  parameters  in  a preprogrammed 
way.  Gain  scheduling  is  easy  to  implement  in  computer 
controlled  systems  provided  that  there  is  support  in  the 
available  software.  This  is  an  ad  hoc  practice  guided  by 
heuristics  rule  of  thumb  [lev], 

A main  problem  in  the  design  of  systems  with  gain 
scheduling  is  to  find  suitable  scheduling  variables.  This 
normally  done  based  on  knowledge  of  the  physics  of  a 
system.  The  concept  of  gain  scheduling  originated  in 
connection  with  the  development  of  flight  control  systems 
(FCS).  In  FCS  applications,  generally  the  Mach  number 
and  the  dynamic  pressure  are  measured  by  air  data  sensors 
and  used  as  the  scheduling  variables.  In  this  paper  only  the 
dynamic  pressure  is  used  as  the  scheduling  variable  since 
the  air  vehicle  to  be  controlled  has  a relatively  narrow 
speed  envelope.  Figure  1.  shows  the  gain  scheduling 
scheme,  which  can  be  viewed  as  a system  with  feedback 
gains  adjusted  using  feedforward  compensation. 


Figure  1.  Gain  scheduling 


Paper  presented  at  the  RTO  SCI  Symposium  on  “Warfare  Automation:  Procedures  and  Techniques  for 
Unmanned  Vehicles”,  held  in  Ankara,  Turkey,  26-28  April  1999  and  published  in  RTO  MP-44. 
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4 Direct  Eigenspace  Assignment 

4.1  Background 

Eigenstructure  assignment  is  a useful  tool  that  allows  the 
designer  satisfy  damping,  settling  time,  and  mode  decoupling 
specifications  by  choosing  eigenvalues  and  eigenvectors. 
Moore  [moo]  in  1976  and  Shrinathkumar  [shri]  in  1978  made 
the  first  discussions  on  this  approach.  Andry  et  al.  [andr] 
have  applied  eigenstructure  assignment  to  designing  a 
stability  augmentation  system  for  lateral  dynamics  of  L-101 1 
aircraft  in  1983.  Sobel  and  Shapiro  [sobel]  used  this  method 
to  design  dynamic  compensators  for  the  same  aircraft  in 
1984.  Later  Sobel  et  al.  [sobe2]  proposed  an  systematic 
method  for  choosing  elements  of  the  feedback  gain  matrix 
which  can  be  suppressed  to  zero  with  minimal  effect  for  and 
F-18  HARV  on  assignment.  Kautsky  et  al.  [kaut]  proposed 
some  robustness  measures  for  system  eigenvalues  to  be  least 
sensitive  to  parameter  variations.  Direct  Eigenspace 
Assignment  (DEA)  method  was  developed  in  1986  by 
Davidson  and  was  used  to  design  lateral-directional  control 
laws  for  F-18  FIARV  of  NASA  in  1992.  This  method  allows 
designers  to  shape  the  closed  loop  response  by  direct  choice 
of  both  desired  eigenvalues  and  eigenvectors.  It  is  called 
“direct”  because  unlike  some  eigenstructure  assignment 
algorithms  feedback  gains  are  determined  in  a single 
iteration.  Davidson  et  al.  [davil]  has  shown  that  during  this 
effort  DEA  has  been  demonstrated  to  be  a useful  technique 
for  aircraft  control  law  design  and  issued  some  guidelines  for 
lateral-directional  flying  qualities  for  high  performance 
aircraft  [davi2]  in  1996. 

4.2  DEA  Formulation 

Given  the  observable,  controllable  system 

x = Ax  + Bu , x e Rn , u e Rm  (4.  i ) 

and 

y = Cx  + Du , ye  R^  (4.2) 

The  total  control  input  is  the  sum  of  the  augmentation  input 
wt.  and  pilot’s  input  ur. 

u = tip+uc  (4.3) 

The  measurement  feedback  control  law  is 

u c = Gy  (4.4) 

Solving  for  u as  a function  of  the  system  states  and  pilot’s 
input  yields 

u = [/,„  - Gd\XGCx  + [lm  -GD\A  up  (4.5) 

The  system  augmented  with  the  control  law  is  given  by 
x = (A  + B[lm  - Gd\ ' GC)x  + B[l m -GD\\ip  (4.6) 


The  spectral  decomposition  of  the  closed  loop  system  is 
given  by 

(A  + B[lm-GDy]GC)vi=livj  (4.7) 

for  /=/ n where  Ai  is  the  i'1'  system  eigenvalue  and  V,  is 

the  associated  i,h  system  eigenvector.  Let  Wj  be  defined  by, 

wi  = [/»,  - GC>\ ' GC\'j  (4.8) 

Substituting  this  result  into  equation  4.7  and  solving  for 
V,one  obtains 

V,-  = [y«  ~AX'Bwi  (4-9) 

This  equation  describes  the  achievable  i'h  eigenvector  of  the 
closed  loop  system  as  function  of  the  eigenvalue  2,  and  the 
eigenvector  tr, . By  examining  this  equation  one  can  see 
that  the  number  of  control  variables  (m)  determines  the 
dimension  of  the  subspacc  in  which  the  achievable 
eigenvectors  must  reside. 

Values  of  ir,  that  yield  an  achievable  eigenspace  that  is  as 
close  as  possible  in  a least  squares  sense  to  a desired 
eigenspace  can  be  determined  by  defining  a cost  function 
values  associated  with  the  i'h  mode  of  the  system 

J,  = t; - v,/(  ]"  £?,/,  [v„  - vlt  ] (4.10) 

for  /=/,...,«  where  v,  is  the  i"  achievable  eigenvector 
associated  with  eigenvalue  2,  , v(/  is  the  i"1  desired 

eigenvector  and  Qt  is  an  n-by-n  symmetric  positive  semi- 

definite  weighing  matrix  on  eigenvalue  elements,  " denotes 
the  complex  conjugate  transpose  operator.  This  cost 
function  represents  the  error  between  the  achievable 
eigenvector  and  the  desired  eigenvector  weighed  by  the 

matrix  Q<t  . 

Values  of  w,  that  minimize  ./,  are  determined  by 

substituting  (4.9)  into  cost  function  for  Vu  , taking  the 

gradient  of./,  with  respect  to  u>  , setting  this  result  equal  to 
zero,  and  solving  for  uy  This  yields 

«’,=k " Q^A.,]' A,,"  Q,,vtli  (4.11) 

where 

A<i,  =[v„-4]  '« 


(4.12) 
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By  concatenating  the  individual  w,’s  column-wise  to  form  W 
and  Va  column-wise  to  form  Va  equation  can  be  expressed 
in  matrix  form  by 

W = [l,„-GDl'GCVa  (4.13) 

The  feedback  gain  matrix  that  yields  the  desired  closed  loop 
eigenvvalues  and  achievable  eigenvectors  is  given  by 

G = w[CVa  + DW\'  (4.14) 

4.3  DEA  Design  Algorithm 

A feedback  gain  matrix  that  yields  a desired  closed  loop 
eigenspace  is  determined  in  the  following  way  [davil]: 


vary  non-linearly  and  tend  to  decrease  in  absolute  value 
with  increasing  dynamic  pressure. 

The  PI  controller  gains  KP1,  Kn,  KP2,  K]2  were  then 
calculated  by  the  Ziegler-Nichols  closed-loop  method  [astr] 
and  lastly,  the  rate  of  climb  adjust  KALT  , determined  by 
trial  and  error  for  the  best  rate  of  climb. 

5.2  Inner  Loop  Controller 

The  microcontroller  based  controller  consists  of  gain 
scheduler,  which  infact  resides  as  a module  inside  the  flight 
management  software,  dynamically  calculates  the  elevator 
demand  in  the  following  scheme: 

1 .  Calculate  air  density:  p = p0(l  -0,00002256h)4256 


1.  Select  desired  eigenvalues  Acl  , desired  eigenvectors 
vd  and  desired  eigenvector  weighing  matrices  Qtl  . 

2.  Calculate  w,’s  using  equation  4.1 1 and  concatenate  these 
column-wise  to  form  W. 

3.  Calculate  achievable  eigenvectors  Va  ’s  using  equation 
4.9  and  concatenate  these  column-wise  to  form  Va. 

4.  Calculate  feedback  gain  matrix  G by  equation  4. 14. 

5 System  Design 


2.  Calculate  dynamic  pressure:  Q = — pV2 


3.  Using  the  gain  equations  that  fit  the  gain  plots  (Figures 
3. a thru  3.d),  calculate  the  gains  Kv  , Ka , Kq  and  K0. 

4.  Calculate  Elevator  demand: 


8e-demand  = 0-demand  + Kv.V  + Ka.a  + Kq.q  + Ke.0 

The  similar  effort  is  performed  for  the  throttle  control  loop. 
Altitude  and  speed  gains  are  calculated,  by  the  controller 
according  to  the  maximum  and  minimum  rates  of  climb  and 
forward  acceleration  specifications  respectively. 


5.1  Methodology 

In  Figure  4.  system  design  is  depicted.  System  states  for 
longitudinal  control  system  are  X = [V,  a,  q,  0,  h];  airspeed, 
angle  of  attack,  pitch  rate,  pitch  angle  and  altitude.  Objective 
is  to  design  both  the  altitude-hold  and  the  speed-hold 
autopilots  for  this  system.  The  related  plant  inputs  are  U = 
[5e,  8t];  elevator  and  throttle  deflection.  Piloted  command 
inputs  are  desired  speed  and  desired  altitude.  Thus, 
theoretically,  four  gains  are  required  for  the  elevator  inner 
loop  feedback,  another  four  is  required  for  the  throttle  inner 
loop  feedback. 

After  determining  a particular  flight  condition,  the  desired 
inner  loop  closed  loop  eigenvalues  and  vectors;  the  DEA 
algorithm  of  section  4.3  was  utilized  in  order  to  calculate  the 
feedback  gains  required  for  the  inner  loop  compensators.  The 
DEA  algorithm  in  section  4.3  was  coded  in  XMATH™  script 
to  enable  the  designer  to  repeat  this  process  for  all  27  flight 
conditions,  within  a few  seconds  for  each  condition.  In  this 
application  dynamic  pressure  Q was  selected  as  the 
scheduling  variable,  which  is  0.5pV2  , where  p is  the  air 
density  in  kg/nr’  and  V is  the  total  airspeed  in  m/sec.  This 
automatic  gain  calculation  process  was  succeeded  by  plotting 
each  gain  versus  the  dynamic  pressure  of  the  corresponding 
flight  condition. 

All  eight  gains  were  plotted  against  dynamic  pressure  and 
these  were  fit  to  polynomials  to  yield  gain  formulas  as  a 
function  of  dynamic  pressure.  Four  of  these  gains  are  shown 
in  Figure  3-a  to  3-d.  Plots  show  that  all  the  longitudinal  gains 


5.3  HIL  Simulation 

Hardware  in  the  loop  simulations  were  performed  on  a non- 
linear aircraft  model  that  covers  the  complete  flight 
envelope  in  MATRIX-X/AC-100™  environment.  Control 
system  of  Figure  4.  was  employed.  Sensor  were  modeled 
and  the  control  system  (the  DSP  processor)  is  added  into 
the  system  loop.  A 100  m climb  demand  was  given  to  the 
system  while  the  vehicle  was  in  equilibrium,  (cruising  at  60 
m/s,  with  trimmed  angle  of  attack  and  pitch  angle  of  1 .8°,  at 
1000  m altitude)  at  1 0th  second.  Altitude,  pitch  angle,  angle 
of  attack,  speed,  throttle  and  elevator  demands  were 
observed.  For  longitudinal  control  design  an  acceptable 
climb  performance  (approx.  4 m/s)  is  achieved  in  the  cruise 
speed  of  60  m/s.  Control  surfaces  move  in  a moderate 
speed  and  within  limits,  (see  Figures  2-a,  2-b,  2-c). 

6 Conclusion 

Direct  eigenspace  assignment  is  used  guarantee  the  inner 
loop  stability  of  the  closed  loop  system.  As  for  longitudinal 
control,  short  period  and  phugoid  modes  are  damped  to  the 
same  nominal  prescribed  frequency  and  damping  factor 
values  for  all  flight  conditions  so  that  compensated  system 
shall  appear  identical  poles  in  all  conditions.  This  gives  the 
designer  the  opportunity  to  design  the  outer  loop  easier 
without  employing  any  further  gain  schedule  for  the  rest  of 
the  system.  Design  philosophy  is  simply,  parametrization  of 
the  inner  loop  gains  via  dynamic  pressure  so  that  the  the 
inner  loop  dynamics  always  appear  almost  the  same  to  the 
outer  guidance  controller  dynamics. 
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Gain  scheduling  had  a positive  effect  on  control  surface 
effectiveness  and  robustness.  DEA  had  positive  effects  on 
overall  stability  of  the  system.  These  two  methods  have  been 
merged  successfully  to  yield  a high  performance  controller 
that  employs  the  advantage  of  each  method. 
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Figure  2.  Controller  in-the-Ioop  system  response  to  1000  m 
climb  demand 


(both  altitude-hold  and  speed-hold  autopilots  ON) 
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Figure  3-a.  Gain  vs.  dynamic  pressure  plot  and  approximate 
curve  fit  with  its  equation  for  Kv 


KA  = 1 E-12Q3  - 3E-08Q2  + 0.0003Q  - 1.1757 


Q 


Figure  3-b.  Gain  vs.  dynamic  pressure  plot  and  approximate 
curve  fit  with  its  equation  for  Ka 


Figure  3-c.  Gain  vs.  dynamic  pressure  plot  and  approximate 
curve  fit  with  its  equation  for  K, 


KTH  = -2E-13Q3  + 3E-09Q2  - 2E-05Q  + 0.0246 


Figure  3-d.  Gain  vs.  dynamic  pressure  plot  and  approximate 
curve  fit  with  its  equation  for  Ke 


Filter 


Figure  4.  Longitudinal  Control  System  Desii 


